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> ■ ABSTRACT 

I Using simultaneous observations in X-rays and optical, we have performed a homo- 

■ geneous analysis of the cross-correlation behaviours of four X-ray binaries: SWIFT 

■ J1753.5-0127, GX 339-4, Sco X-1, and Cyg X-2. With high time-resolution observa- 
I tions using ULTRACAM and RXTE, we concentrate on the short time-scale, St < 20 s, 

00 ' variability in these sources. Here we present our database of observations, with three 

' simultaneous energy bands in both the optical and the X-ray, and multiple epochs of 

, observation for each source, all with '^second or better time resolution. For the first 

time, we include a dynamical cross-correlation analysis, i.e., an investigation of how 
the cross-correlation function changes within an observation. We describe a number 
of trends which emerge. We include the full dataset of results, and pick a few striking 
' relationships from among them for further discussion. 

I We find, that the surprising form of X-ray/optical cross-correlation functions, a 

positive correlation signal preceded by an anti-correlation signal, is seen in all the 
sources at least some of the time. Such behaviour suggests a mechanism other than 
reprocessing as being the dominant driver of the short-term variability in the optical 
emission. This behaviour appears more pronounced when the X-ray spectrum is hard. 
Furthermore, we find that the cross-correlation relationships themselves are not stable 
in time, but vary significantly in strength and form. This all hints at dynamic interac- 
tions between the emitting components which could be modelled through non-linear 
or differential relationships. 

Key words: Stars: neutron; Stars: black holes; X-rays: binaries; 



1 INTRODUCTION 

X-ray binaries (XRBs) are systems in which an ordinary star 
and a compact object are in close orbit. Matter falls onto the 
compact object and is accreted by it. The in- falling material 
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releases its binding energy, powering energetic radiation and 
matter out-flows. In the most typical scenario (and the only 
one considered here), the donor star is an evolved star at 
the end of its main sequence life, and over-filling its Roche 
Lobe. The compact object can be a black hole or a neu- 
tron star; white dwarf systems are not normally considered 
XRBs, since the majority of the luminosity is released at 
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somewhat lower energies, corresponding to the larger radius 
at which in-falling material is halted. See lLewiri et all (|l997l ) 
for a description of X-ray binary phenomenology and theory, 
and possible models for the accretion process. 

Since XRBs are most manifest in the X-ray band, much 
observation of them has focused here. This is particularly 
true of timing investigation, since X-ray instrumentation has 
been capable of high temporal resolution since the start. A 
wide range of time-domain behaviour have been catalogued 
for the various binary systems, such as broad-band timing 
noi se and a v ariety of quasi-periodic oscillations (QPOs). 
See Ivan derl^ lis ( 2004) for a review of X-ray timing obser- 
vations. 

The X-ray emission corresponds to the bulk of the lu- 
minosity and should come from the smallest radii (and cor- 
respondingly highest particle energies) around the compact 
object, where the dynamical time-scale is the shortest. This 
strongly suggests, therefore, that any optical variability on 
a similar time-scale to the X-rays must occur by reprocess- 
ing, i.e., the absorption and re-emission of X-ray radiation 
by cooler material at greater radii in the accretion disc or on 
the surface of the donor star . Alth ough this idea is not new, 
Ivan Paradiis fc McClintockl (|l994) modelled this in detail, 
and presented evidence that the optical luminosities were of 
a similar magnitude to what might be expected, if repro- 
cessing were dominant in the optical. 

Optical emission by reprocessing of higher-energy radia- 
tion, e.g. X-rays, has undoub t edly been seen in X-r ay bursts 
(see e.g.. lGrindlev et al ]|l978l : iPedersen et "allll982l ) . In these 
cases, it is believed that the burst originates in an explosive 
event on, or immediately above, the surface of the neutron 
star. The optical emission is seen to rise sharply very soon 
(<ls) after the X-ray burst, and can be well-modelled as 
absorption and re-emission of the original flare, mainly in 
the accretion dis c at larger radii, a nd also from the donor 
star surface (e.g. iHvnes et al] l|2006l )). The cross-correlation 
function (CCF) for this case is a single positive peak centred 
at positive lags. The X-ray light-curve in such a scenario re- 
flects primarily the temperature of the neutron star surface 
and its short-term behaviour is not connected with the ac- 
cretion disc or jets. Nevertheless, X-ray bursts serve to show 
that the optical emission, and its rapid variability, can be 
dominated by reprocessing at least for a short time, given 
certain specific circumstances. 

Reprocessing has also been demonstrated through the 
technique of Doppler Tomography for individual spectral 
lines. The known period of the binary orbit is used to find 
the phase and velocity amplitude of emission components 
within a strong (emission) line. This can be transformed 
to map the emission lo cations within the binary system 
(see[C asares et al.ll200^ . for an example of this applied to 
a low mass X-ray binary). Similarly, echo tomography has 
been attempted in a number of cases using the fast vari- 
abilit y near the co mpact object and the echoes this can in- 
duce (|Hvnesll2005f ) . especially for larger-scale accreting sys- 
tems , the active galactic nuclei (AGNs, e.g., iRaiteri et al.l 
I2OO8I ). In contrast to the optical, infra-red emission may 
well be contaminated or dominated by a jet: a relativis- 
tic, energetic outflow originating near the compact object 
and travelling along the polar axes. The emission from a 
jet is therefore causally connected with the central accre- 
tion engine, but fluctuations should happen later than, and 



in response to, the higher-energy em ission in the core. See 
for example lEikenberrv et al. ll2008h for the case of a galac- 
tic micro-quasar, or lCasella et al.l 1 2010l ) for the most recent 
such observations of one of the systems under consideration 
in this paper. 

Optical high time resolution observation is still a rela- 
tively young field, since CCDs capable of fast readout with 
low noise are relatively new, and the previous generation of 
detectors suffered from low quantum efficiency and/or high 
dead-time. Due to the difficulty of coordinating the schedul- 
ing of multiple observatories, not many simultaneous sub- 
second time-resolution, multi-band observations exist in the 
liter ature. 

iKanbach et al.l (|200ll ) observed the galactic black hole 
XRB XTE J1118-I-480 simultaneously with the Rossi X- 
ray Timing Explorer (RXTE) and OPTIMA, an optical 
single-pixel fast photometer, observing white light. They 
found two surprising results in the correlation analysis: 
the optical auto-correlation function (ACF) had a smaller 
width than X-rays, and the CCF had a small component 
which was negative and with negative lag (i.e., optical ar- 
riving before X-rays). Neither of these factors would be ex- 
pect ed for the simplistic r eproc essing model. In further anal- 
ysis, ISpruit fc KanbachI (|2002h commented that the cross- 
correlation function was not stable through the observation, 
and contained two statistically separable principal compo- 
nents, both of similar dip-and-peak shape, but differing in 
time-scale. A simplistic explanation was presented to explain 
the precognition dip, of a local density and vertical height en- 
hancement in the disc obscuring the inner X-ray region, and 
then causing increased X-ray emission some seconds later 
as the blob is accreted. No numerical or detailed modelling 
was done, however. Further discussion on the possible phys- 
ical processes and time-scales in the system was provided by 
Hyncs ot al^ (2003)^ 

iDurant et all lj2008l ) and lCandhi et al.l ('2OO8) increased 
the number of fast {5t <ls) optical/X-ray cross-correlation 
analyses in the literature to three. They also suggested that 
the fast functions seen were not compatible with reprocess- 
ing as the dominant source of optical variability. That data 
is included in our analysis below. 

In summary, the reprocessing model is able to explain 
the energetics of the optical emission for X-ray binaries in 
general, and many of the details of the spectrum and long- 
term behaviour that would be expected. It has problems, 
however, with the fastest time-scale variations, and it is 
these that we wish to concentrate on in this paper. Time- 
scales < 2s are shorter than the light-crossing time for a 
typical accreting binary system (depending, of course, on 
the binary separation). 

Here we analyse multiple ULTRACAM/RXTE obser- 
vations for four XRBs: SWIFT J1753.5-0127, GX 339-4, 
Sco X-1, and Cyg X-2. For each observation, we consider 
three energy bands in the optical and three in the X-rays, in 
a homogeneous way across the dataset. Also, we investigate 
the dynamical behaviour of the CCF, a technique which is 
not commonly considered. The total number of presentable 
relationships are large; we include the full set in Appendix 
A, and we draw out a few examples to illustrate the types of 
behaviour we find in the data. What we see suggest that a 
large fraction of the fast variability seen in optical emission 
is not produced by reprocessing. Therefore, some explana- 
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tion is required for why reprocessed variability is not more 
pronounced. 



2 OBSERVATIONS 

Table [T] lists the observations considered as part of this si- 
multaneous observation compilation. As shown, some of this 
data has been presented before in some form. 

ULTRACAM is a high-speed optical camera, which in- 
corporates dichroics to split the input beam into three en- 
ergy bands. The detectors are fast read-out CCDs with read- 
out noise of only 3.5 e—. Masked data buffer areas on each 
chip enable an image to be very quickly shuffled from the 
active region and then to be read while the next is expos- 
ing, reducing the dead-time to near zero. By reading only 
small windows around the target and reference star, sam- 
pling up to 500 Hz is possible. A dedicated GPS receiver 
provides reliable time-stampin g to sub ms acc uracy. See 
iDhillon fc MarshI |2001) and D hillon et al.1 (120071 '). We anal- 
ysed each dataset using the custom ULTRACAM pipelin^, 
using variable-aperture photometry relative to a bright field 
star. Slow variations due to the difference in colour of the two 
stars make negligible effect on the time-scales under investi- 
gation. For these observations, ULTRACAM was mounted 
on either the 8.2 m Very Large Telescope (VLT), Paranal, 
Chile, or the 4.2 m William HerscheU Telescope (WHT), La 
Palma, Spain. 

The dichroics in ULTRACAM split the light into 
the blue, middle and red parts of the spectrum, but fil- 
ters are typically used before each detector. For SWIFT 
J1753. 5-0127 and GX 339-4, broad r'g'u' fihers were used, 
but for Sco X-1 and Cyg X-2, narrow filters were used 
in the r-continuum (i.e., a spectral region without lines), 
and in the vicinity of the Bowen/He-11 blend (the broad u' 
was retained). The purpose was to try to isolate the emis- 
sion in the narrow emission lines, but the majority of the 
flux in each filter is still dominated by the continuum (see 
iMunoz-Darias et al.ll2007l ). 

RXTE is an X-ray observatory optimised to tim- 
ing work. The Proportional Counting Array (PCA, 
Ijahoda et al.|[l99^ ) instrument has a large effective area and 
is sensitive to photons with energies 2-60 keV. Here we only 
consider the more sensitive 2-20 keV range. Although a num- 
ber of data collecting settings are possible, the intrinsic time 
resolution is of the order /is with full energy resolution for all 
but the brightest sources. The limiting factor is the teleme- 
try bandwidth. The data were processed using the standard 
FTOOLS package^, and three background-subtracted light- 
curves produced for each dataset: 2-4, 4-8, and 8-20 keV. 
Note that for a non-imaging instrument such as the PCA, 
background subtraction is in general a complicated issue; 
estimates can be based on empirical measurements, and de- 
pend on the spacecraft attitude and recent trajectory. We 
utilised the latest South Atlantic Anomaly passage history 
and background models, 2009. For such bright sources how- 
ever, these make little difference, and the noise at time-scales 
or order 1 s or less is dominated by source photon statistics. 



3 ANALYSIS 

Our aim in this paper is to compare the cross-correlation 
functions between the observed X-ray and optical flux for 
each source. This is a simplistic determination of the rel- 
ative time- variability of the two wavebands, i.e., to what 
degree do they vary in similar ways, and what is the relative 
timing of these variations. It provides a clear, uncomplicated 
comparison of two signals with respect to time, but is by ne- 
cessity an averaging process. To go from a significant signal 
in a CCF to a realistic picture of what the two light-curves 
are doing is not always straightforward. We endeavour to de- 
scribe the typical behaviour that gives rise to the CCFs seen. 
We note that alternative methods for time-series compari- 
son exist, such as the coherence and Fourier cross-spectrum, 
but in our experience these are no more transparent, and 
indeed are intimately related to the CCF. We only consider 
CCFs from here on. 

We calculated CCFs for each simultaneous observation 
by first re-binning the light-curves to the best time reso- 
lution of the poorer of the pair. CCFs were calculated in 
windows of size 50 s and stepping 25 s, i.e., overlapping sec- 
tions where every second CCF in independent. A mean CCF 
is then derived for each time-series pair for each observation. 
This acts effectively as a high-pass filter, and variations with 
time-scales longer than ~50s will not appear. For the pur- 
poses of this paper, we are only interested in the fastest 
variability, since this must be predominantly due to fluctua- 
tions arising close to the compact object, where the dynam- 
ical time is the shortest. 

All the CCFs were considered, for each object, obser- 
vation and pairs of energy bands. In Section U we present 
a selection of average CCFs, pointing out some similarities 
between the binary systems, and some trends in the data. 
The full set of CCFs is presented in Appendix A. Note that, 
although a comparison star in the field is used in each ob- 
servation to mitigate atmospheric effects, the weather still 
greatly a ffects the signal qua lity achieved in a CCF, see for 
example (|Gandhi et al.ll2010h and the comparison of succes- 
sive nights, where we find that the CCF features are most 
pronounced on the night with best weather 2007-06-18. 
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3.1 Dynamic analysis 

In the initial analysis above, the CCFs were the average 
of many 50 s sections throughout an individual observation. 
These segments can, however, be interesting independently. 
A convenient display for such a set of CCFs, is to stack them 
above each other, with one line representing a single CCF, 
and colours along that line giving the strength of the correla- 
tion for a given lag, at a given point in the observation. Such 
a representation is similar in style to the more common dy- 
namic Fourier-gram. Note that because the segments over- 
lap, adjacent CCFs are not independent, but every second 
one is. This produces slightly smoother-looking output, and 
is intended to catch strongly correlated events that would 
am/JtijiSHiBfeiaiitaB the border between two segments. No ad- 
ditional numerical smoothing has been applied. 
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Table 1. Log of simultaneous ULTRACAM/RXTE observations 



Object 


Date 


ULTRACAM 


Optical time 


Simultaneous 


Reference^ 






site 


resolution (s) 


coverage (s) 




oWir 1 JliOo.o — UiZi 


ZUU (-UD-1Z 


V ijl 


U.i4 




JJurant et al. (^zUUo } 




2008-08-20 


WHT 


0.2 


2400 




GX 339-4 


2007-06-14 


VLT 


0.14 


2800 






2007-06-16 


VLT 


0.13 


2700 


Gandhi et al. ("20081 




2007-06-18 


VLT 


0.05 


3300 


Gandhi et al. ("20081 


Sco X-1 


2004-05-17 


WHT 


0.1 


950,760,730 






2004-05-18 


WHT 


0.25-0.5 


260,700,950 


Muiioz-Darias et al. ("2007) 




2004-05-19 


WHT 


0.3 


115 




Cyg X-2 


2007-10-16 


WHT 


2.0 


880 






2007-10-17 


WHT 


2.02 








2007-10-18 


WHT 


2.02 








2007-10-19 


WHT 


2.0 


600,710 






2007-10-20 


WHT 


2.0 


530 






2007-10-21 


WHT 


2.0 


1050 





^ where no reference is given, the data are presented here for the first time. 
^ the data quality in these observations is too poor for consideration. 



3.2 CCF noise 

Aside from intrinsic changes in the signal, two types of noise 
also contribute to the differences between any two CCFs and 
to the structures within a single CCF. The first is simple: un- 
correlated (white) noise inherent in any measurement. This 
is dominated in our case by the Poissonian statistics of pho- 
ton counting, but read noise can also contribute. Since the 
cross-correlation calculation is a simple multiply-and-sum 
process, these errors are easy to propagate. Longer time- 
scale noise is also present in the data, however. These are 
due both to slow systematics (e.g., seeing variations) and 
red noise processes intrinsic to the emission processes them- 
selves, and act on a range of time-scales. These make the 
estimation of the noise contribution to CCFs rather more 
complicated. 

One simplistic method to establish the significance of 
an average CCF signal might be to measure the standard 
deviation of all values for a given lag, compared to the aver- 
age value. This would only be valid in the case that the CCF 
is assumed to be stationary, but in our dynamic analysis we 
find this generally not to be the case, so we do not consider 
this method further. 

The first-order approximation we apply, to see whether 
there is any significant signal in a given CCF, is to empiri- 
cally measure the noise in the CCF itself, given by the me- 
dian of the difference between successive points. This value 
is a combination of various noise sources, and as such, a 
useful estimator for analysing a given CCF. Signals falling 
more than five times the value are likely to be significant. 
In the averaged CCFs shown in Appendix A, we give both 
a raw CCF scale, and a a scale (CCF/a), so that the y-axis 
becomes a rough estimate of significance. For any significant 
average signal, one can inspect the corresponding dynamic 
CCF plot, to see where in the observation the signal mostly 
originates from; it may be steady throughout the observa- 
tion, intermittent, or changeable in height and shape. 

Having established a candidate CCF signal, one can ask 



how well the features within it can be localised, This we ap- 
ply to the dynamic CCFs, since we find the possible variation 
in peak height and location interesting. 

In order to estimate the confidence in establishing the 
amplitude and position of a particular CCF peak, we follow 
the empirical estimation of ISmith fc VaughanI (|2007l ) who 
sh ow a simple imp lementation of the theoretical derivation 
bv iBartlettI ()l955h . This method considers both the corre- 
lated and un-correlated noise in each light-curve, as empiri- 
cally sampled by the auto-correlation functions. It provides 
the uncertainty on each point in a calculated CCF. Since 
the method is only applicable to a stable CCF relationship, 
our approach here is necessarily simplistic. To estimate the 
confidence interval on an individual peak, one can fit a sim- 
ple function (e.g., a Gaussian) with a standard chi-squared 
process incorporating these uncertainties. Such a fit is only 
a good estimator for the case that one has good reason to 
believe the Gaussian function fit is a good description of 
the data. The CCF functions are clearly more complicated 
than this, but we nevertheless think this approach is a good 
approximation for characterising peaks. 



4 RESULTS 

We include the complete set of cross correlations, averaged 
and dynamic, in Figures Al to A34, in Appendix [X] They 
serve to illustrate, that the range of behaviours that we find 
are not particular to a small subset of the sample, but gen- 
eral features. The figures are shown with consistent colour 
mapping. For each averaged CCF, we give both significance 
levels, based on our empirical estimate (see above), and ab- 
solute CCF values. In this Section, we focus on a few inter- 
esting CCFs, in order to best illustrate some of the trends 
in the data. 

In Figures[T]to|3]we present a selection of the CCFs cal- 
culated for our sample of objects, energies and epochs. The 
style of observation and the analysis have been conducted 
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in a homogeneous manner. We extract a few trends that are 
apparent in the data. First we present typical CCFs for each 
source, then look at the trends with X-ray and optical en- 
ergies, and finally display some examples of the dynamical 
analysis. The latter requires a little more careful consider- 
ation, since the changes seen in the CCF through a single 
observation could in principle be caused by noise processes 
alone - see Section 14.31 Although in Section [5] we attempt 
to draw conclusions from what we find, and speculate on 
the physical scenario behind them, principally the results 
below stand as a phenomenological study, for interpretation 
in future theoretical works. 

For the purposes of these example CCFs, we have picked 
plots that show particular features clearly. The full set are 
included in Appendix A, and the reader is encouraged to 
look through these to see the frequency with which the var- 
ious features and trends outlined below occur. 

4.1 Anti-correlations 

Figure[T]shows example CCFs for each of the four binaries of 
our sample. They show the following typical form: negative 
correlation at negative lags {CCF < 0, St < 0) and some- 
what weaker positive correlation at positive lags. The sharp 
slopes seen around St = indicate that in each case, much 
of the signal is due to variations at times-scales comparable 
to the the shortest measured. Even shorter time-resolution 
might of course be present, but is unobservable in these data. 
The most important point, is that of the four systems ob- 
served, all showed a similar CCF form at least some of the 
time. 

Not all of the observations show the typical CCF form 
described above. In a number of observations, the level of 
noise, predominantly Poissonian (with some contribution 
from, for instance, seeing and transparency variation), make 
any CCF signal relatively weak or unobservable. CCFs show- 
ing a single, strong positive peak at positive lags, as would 
be expected from reprocessing, are notable by their rarity. 
Furthermore the optical ACFs (not shown) are consistently 
of similar width or sometimes narrower compared to the X- 
ray ACFs. It would thus seem that the majority of the fast 
variability in the optical is not caused by reprocessing. The 
existence of a CCF signal indicates, however, that the two 
lig ht-curves are n ot independent, but connected. As noted 
bv ICandhi et al.] (|2008), it may be possible to describe a 
complicated dip-and-peak CCF as a broader trough and sim- 
ple peak, both centred near St = 0. This was motivated for 
the case of GX 339—4 by the similarity between the width 
one would ascribe to the broad trough and the width of the 
optical auto-correlation function. 

The sources presented here have various emission states 
in their epoch to epoch evolution, typically classified by their 
location on an X-ray flux versus hardness (i.e., flux ratio 
of a higher and lower energy band) diagram. The obser- 
vations here all refer to times when the sources were in a 
low-hard state or flaring states, with a luminosity well be- 
low the maximum known, and the spectral slope in the X- 
rays harder. For example, SWIFT J1753.5-0127 has since 
its discovery, never been seen to leave the low-hard state 
l|Durant et al.ll2009l '). For Sco X-1, where the hardness of the 
source changed markedly during the observation campaign, 
the strongest CCF signatures correspond to times when the 



X-ray emission was hardest, although this generally also co- 
incided with times when the overall flux was the highest. 
Note that none of the systems was seen in t he classical 
High /Soft state during these observations (e.g.. Id'Ai et al.l 
I2OO7I for a typical evolution between various states for Sco 
X-1); in the high/soft state, one expects copious thermal 
emission to dominate the lu minosity, and fast vari ability to 
be much less apparent (e.g.. iBradshaw et al.|[2007l ). 



4.2 Energy dependence 

Figure [2] shows two examples of a systematic trend in 
the CCFs with X-ray energy band. For the Cyg X-2 data 
(left panel), one sees that the positive component becomes 
stronger and peaks closer to zero lag for successively higher 
energies. The change in strength, from ~0 to 0.2 is dramatic. 
These curves were generated from the same observation. The 
right-hand panel for Sco X-1 shows a similar result: the curve 
becomes stronger and peaks closer to zero lag with higher 
energy, while the negative component shows little change. 
The change in correlation strength in this case is more sub- 
tle. 

Figure [3] shows two examples of a systematic trend in 
the CCFs with optical band. For GX 339-4 (left panel), 
one sees that both r' and g' show the same structure, anti- 
correlation near zero lag, except for a sharp peak very close 
to zero lag, but that the relationship is stronger for r', par- 
ticularly for the positive part, u' shows something different, 
broad negative correlation at negative lags and positive cor- 
relation at positive lags. Note that, due to the lower time- 
resolution of the u' data, the sharp peak a,t St — for r' and 
g' is not excluded in the u' CCF. The right-hand panel for 
SWIFT J1753. 5-0127 shows a different result: the g' CCF 
curve has the same shape, but stronger than the r' one. 



4.3 Dynamic behaviour 

Figure U shows three examples of the results of our analysis. 
Each figure represents a set of CCFs from two time-series in 
one particular observation. Each horizontal line is a single 
CCF as measured at the time T in the observation. This time 
is given as an offset since the start of simultaneous coverage. 
The colours along the horizontal line represent the strength 
of correlation for different values of lag, St, along the x-axis. 
The colour scales are given for each panel individually, and 
can be compared with the values on the y-axes of Figures [1] 
toO which are average CCFs, or to the whole set of CCFs 
in Appendix A. 

By eye, a number of features are apparent: 

• for GX 339—4, the positive component seen in the 
average CCF appears for very short bursts and relatively 
strongly, with a low duty cycle through the observation. 
The negative component appears to be present whenever 
the positive component is not. 

• for SWIFT J1753. 5-0127, the negative component is 
ever-present and strong, but the centroid appears to shift, 
and somewhat cross to St > 0. Likewise, the weaker positive 
component is sometimes centred at St < 0. The strength and 
width of both components changes markedly, but without 
any obvious trend through the observation. 



6 Durant et al. 



o 

o 



0.7 



0.6 - 



0.5 - 



0.4 - 



0.3 



0.2 - 



0.1 



0.0 



-0.1 



-0.2 



SWIFT J 1753 
SCO X-1 
Cyg X-2 
GX 339-4 




-20 



-10 





St (s) 



10 



20 



Figure 1. Example average CCFs for our sample of XRBs, demonstrating the ubiquity an anti-correlation component at negative lags 
in each case. These are each in the optical red band and 2-4koV for X-rays. There is a vertical offset of 0.15 between successive curves, 
for clarity. 



• for Sco X-1, the form of the CCF seems to be constant, 
but the strength and location of the peak changes remark- 
ably throughout the observation. 



The dynamic CCF diagrams clearly show a lot of vari- 
ation with time. One must consider whether such variations 
could be caused by the noise processes inherent to the light- 
curves under analysis. Here we present the application of 
such an investigation to a single dynamic CCF diagram 
which shows marked variations, as an illustration. We wish 
to test the significance of the peak variations for case of Sco 
X-1 (2004-05-18), because of all the observations, this is the 
one that shows the most spectacular and unexpected CCF 
variations. Using the Bartlett formalism described in Sec- 
tion 13.21 to assign uncertainties to every CCF bin value, we 
have performed fitting to find the peaks in each line of the 
dynamic CCF. We fitted a single Gaussian function to the 
positive component of each CCF line and minimised the 
statistic. The uncertainty of the location of that peak is then 
the value at which increases by 1. 

The result of such an analysis are shown in Figure [5] 
We find that the shift in the location of the positive peak is 
highly significant, and a peak lag of (5t < at some times is 
marginally significant. Note that here we are assuming that 
a single Gaussian function is a reasonable description of the 



data, which is clearly not the case; it is sufficient, however, 
for illustrative purposes. 

Mufioz-Darias et al.'s (2007) analysis was performed on 
a 180 s window centred at around T — 450 s on this diagram. 
Had we restricted our analysis to just this section of the 
data, we would also have concluded that the major feature 
of the CCF is a positive peak with a centre at around St = 
5. .10 s. Note however that these authors were attempting to 
measure the differences between CCFs, as a way to probe 
the emission lines which are present in the spectral passband 
of one of the filters but not the other. 

Using a similar process for the rest of the dynamical 
CCFs, we find that the features noted above for Figure [3] 
are significant. 



5 CONCLUSIONS/DISCUSSION 

We have presented a number of optical/X-ray CCFs from si- 
multaneous observations of four XRBs, and deeper analysis 
of the energy- and temporal dependence of these. The most 
obvious and important result is, that for such a heteroge- 
neous sample, optical-leading anti-correlations appear to be 
typical of the emission behaviour. Currently the accepted 
standard picture is that the optical emission arising from 
X-ray binaries is dominated by reprocessing of the X-rays in 
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Figure 3. Example average CCFs for single observations (left: GX 339-4, 2007-06-16, 2-4 keV; right; SWIFT J1753.5-0127, 2008-08-20, 
4-8 keV) demonstrating trends with optical band. 



the outer accretion disc, or on the surface of the donor star. 
One would then expect that the cross-correlation of opti- 
cal and X-ray flux shows that any variability in the X-rays 
is mirrored in the optical flux, but arriving to the observer 
a little later to account for the light-travel time across the 
binary system and the absorption/emission time-scale. In 
other words, the CCF should show a single strong, positive 
peak at positive lags {5t > 0, optical arriving after X-rays.). 

Our observations show that the CCFs one typically sees 
are more complicated, and in particular, that signals of anti- 
correlation and optical-leading are ubiquitous, at least in the 
relatively low states in which our binaries were observed. It 
is thus difficult to reconcile the origin of the fastest vari- 
ability in optical being due to X-ray reprocessing. This sug- 
gests that there are additional processes involved in accre- 
tion, that are not part of the standard model, the accretion 
process is more complex than previously thought. 

This does not mean, however, that reprocessing is ab- 
sent altogether: if X-rays are illuminating cooler material. 



it is natural to ask why the reprocessing CCF signal is not 
seen. One simple answer could be that whatever is giving 
rise to the signals we observe is swamping the reprocessing 
signal; but it could also be that the range of light-travel 
times or the reprocessing time-scale is longer than previ- 
ously thought, and that therefore the variable component of 
reprocessing is washed out on the time-scales we are probing. 
Decomposing the CCF into its fourier components an d phase 
lags may be a way to get a round this difficulty ( Malz ac et al.l 
I2OO3I , ICandhi et al.|[201ol '). Certainly, there are strong fluo- 
rescence lines in some of t he spectra, e.g., known since as 
early as iMook et all (ligTll). but ri otably absent for SWIFT 
J1753.5-0127 (|Durant et al.ll2009l '). Using fast spectroscopy, 
it would be possible to measure the variability in lines com- 
pared to X-rays, as opposed to the variability in the contin- 
uum (or the relative variability of lines and continuum). 

Since the column density required to reach an optical 
depth of 1 (i.e., attenuation by a factor e) for photon en- 
ergies >2keV is of the order > lO^^cm"^, and higher for 
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Figure 4. Example dynamic CCFs (upper left: SWIFT J1753. 5-0127, 2007-06-12 (93119-02-02-00); upper right: GX 339-4, 2007-06-18 
(93119-01-04); bottom: Sco X-1, 2004-05-18 (90020-01-02-02); all r/2-4keV, Obs IDs in parentheses). Each horizontal line is a single 
CCF from a 50 s segment, with the colours representing negative (blue) or positive (red) correlation. Successive lines overlap, see text. 



more energetic photons, it could be simply that the pene- 
tration depth is so deep th at the optical re p rocess ed hght 
takes time to emerge (e.g. iMcGowan et~all (|2003[ )1. With 
the X-ray spectra moderately hard as here, a far smaller 
proportion of the luminosity is released in the thermal, 0.5- 
2 keV X-ray band than for soft sources. Conversely, irradia- 
tion by higher-energy photons may be expected to produce 
a str onger temperatu re change in the surface layers of this 
disc (|Wu et al.ll200l] ). It would be interesting to do a sim- 
ilar study at thermal energies (<2keV), and also at much 
higher energies. Furthermore, it is not necessarily justified 
to assume that the X-ray light-curve reaching us, is the same 
as that 'seen' by the material which is producing the optical 
light reaching us. A detailed three-dimensional model would 
be required to describe this process; however we have found 
that the light-curves are related, but not as we might have 
expected. 

In Section 1 we noted that for the particular circum- 
stances of an X-ray burst, optical emission by reprocessing 
is clearly observed. We also mentioned the use of Doppler 
Tomography to map the sites of emission using the varia- 
tions of line profiles with orbital phase, for particular, strong 



emission lines. Such an analysis shows that both some of the 
short-term and of the persistent optical flux originates in 
the accretion disc, in-fall stream, surface of the donor star 
etc. The total amount of flux contained within the emission 
lines is typically small compared to the continuum, and one 
would not expect the line flux to dominate even in a nar- 
row filter centred around a stronger line. Given that there 
are typically not many such fiuorescence lines in the spec- 
trum, one would certainly not expect the flux measured in 
a broad-band filter to be affected by lines. Therefore, the 
fast variability measured here is indicative of the continuum 
emission and not of the fluorescence lines, and there is no 
contradiction between the majority of the CCF signal being 
incompatible with reprocessing, and the lines which are in- 
deed necessarily produced by reprocessing. Furthermore, we 
have not considered slow variability {St > 20 s) here at all, 
where reprocessing is known to occur. 

Ivan Paradiis fc McChntockl (|l994 ) made some theoret- 
ical predictions on the amount of luminosity to be expected 
from reprocessing in an X-ray binary. They propose a linear 
relationship for the visual luminosity Ly oc l}^P^^^ (P is 
the binary period and Lx the X-ray luminosity responsi- 
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Figure 5. Evolution of the positive correlation component of the dynamic CCF for Sco X-1, 2004-05-18 fFigure[4l lower panel). A single, 
positive Gaussian has been fitted to each line, as a way to estimate the centre of the positive correlation peak. The uncertainties on this 
fit have been found using realistic errors (see text) and the error-bars represent l-tr confidence limits. On the left panel, we show the 
best-fit location in &t of the Gaussian, and it's uncertainty, and on the right panel we show the fitted amplitude and its uncertainty for 
each CCF line. Where the marker and associated error bar are missing, the Gaussian fit failed. 



ble for reprocessing). In a linear fit to several well-studied 
LMXBs, they find absolute magnitude 



Mv = 1.57(24) - 2.27(32) log 



where L Edd is the Eddington luminosity. Applying the same 
method, and making reasonable assumptions for the un- 
known Lx/LEdd and using a range of possible distances, 
we have plotted the four objects in our study, plus XTE 
J1118-I-480 next to the relationship given in Equation[T] see 
Figure [5] (first discussed by , Gandhi et aL , 2010, for the spe- 
cific case of GX 339—4). For the purposes of this plot, the 
luminosity is in the 1-100 keV range (as extrapolated from 
spectral fits), and assumed all to contribute to reprocessing, 
and the magnitudes are corrected for reddening, which is an 
additional source of uncertainty. Only for Cyg X-2 would 
one expect reprocessing to be a significant fraction of the 
optical emission. Perhaps it should not be a surprise, then, 
that reprocessing does not appear to make much of a differ- 
ence to the CCFs; but conversely a new source of the opti- 
cal variability, and its connection to X-rays is required. Note 
that the systems showing the most unusual timing properties 
(and the fastest variability), XTE J1118-I-480, GX 339-4 
and SW IFT J1753.5-0127 are also the on es lying furthest 
from the lvan Paradiis fc McClintockl (| 19941 ) relation. It is in- 
teresting to note that these are black hole systems, whereas 
the neutron s tars s ystems lie much closer to the relationship. 

iGandhil (|2009l ) found that a relationship exists between 
the RMS and flux for the three binary systems GX 339—4, 
SWIFT J1753.5-0127 and XTE J1118-h480. The two quan- 
tities can be adequately described by a linear trend (with 
significant scatter), that does not pass through zero, for vari- 
ous time-scales. Such behaviour had been found in X-rays by 



Cyg)« 



Figure 6. Expected a b solute magnitudes predicted by 
Ivan Paradiis fc McClintockl l|l994l ) (solid line) compared to the 
measured magnitudes for the X-ray binaries of interest (mark- 
ers) . The dotted lines indicate a range of plausible distances for 
each object. S is as defined by these authors, the expression in 
square brackets in Equation [T] 



lUttlev et ai]l|2005h . and implies that the light-curves cannot 
be described by a simple shot process of independent fiares, 
but rather that the size and duration of events is related, and 
probably connected to their location in the accretion disc. 
The emission thus cannot originate from many independent 
launching sites. Thus any theory wishing to describe the 
CCF behaviour above must also consider this important re- 
sult. 
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In our observations, variability is seen down to the 
shortest time-scales sampled for every source, and the result- 
ing CCFs have such sharp gradients that some of the corre- 
lation must occur at these time-scales or faster. One wonders 
how much more structure would be discovered with higher 
time resolution, high signal-to-noise observations. Variations 
of the CCFs with X-ray or optical band is an indication 
that the varying component responsible for the CCF sig- 
natures has a different spectrum to the average spectrum. 
Differences between the CCFs of the various bands will be 
affected by the difference signal-to-noise of the individual 
light-curve, differences that are hard to model. 

A number of q ualitative mode ls have appeared since the 
initial discovery in lKanbach et all (|200ll ) of the unexpected 
correlations for XTE J 1118-1-480:- 

• The first att empt was by the same authors in more 
detailed analysis l|Spruit fc Kanbachll2002l ). Their simplistic 
description was of blobs of dense material in the accretion 
disc, having a higher density and therefore larger vertical ex- 
tent. As they travel towards the central source, they would 
partly obscure the X-ray light-source, just before being ac- 
creted and causing an increase in X-ray luminosity. Unfor- 
tunately, the model did not include detailed predictions. We 
would expect, however, either X-ray energy independent be- 
haviour (for partial covering by a dense, opaque blob) or 
variability much stronger in softer X-ray bands (for moder- 
ate optical depth) . 

• iMalzac et all (|2004l ') produced an alternative toy model 
to explain XTE J1118-|-480's CCF. This involves a magnetic 
energy reservoir, where the energy input is a stochastic pro- 
cess, and the output can be channelled either to X-ray lu- 
minosity or the jet. Thus the emission at a moment depends 
not only on the current accretion rate, but on the past his- 
tory; the emission modes are therefore linked in a differential 
relationship. Such a model can quantitatively explain the 
cross- and auto-correlation function, although how to main- 
tain such an energetic magnetised zone is not addressed. 
Coupling between the different emission zones, at various 
time-scales is r equired by the fl ux-RMS relationship men- 
tioned above. (|Yuan et al.l 120051 ) developed a quantitative 
model along the ideas of multiple contributing components 
to t he optical flux. 

• ICasella fc Pe ^ (|2009l) provide yet another explanation, 
this time questioning the assumed synchrotron spectrum ex- 
pected from a magnetised outflow. Typically a broad-band 
flat spectrum is assumed, breaking and falling in the 
IR-range; but this assumes a special magnetic field/outflow 
configuration, developed to match flat radio spectra. For 
suitable magnetic fields, it can be that the synchrotron spec- 
trum had a bump near the optical range, but also provides 
the bulk of the X-ray flux. Thus the emission in the two 
waveband is intrinsically connected, and the relative timing 
depends on the magnetic field strength and geometry in the 
jet. 

• Conversely, part of the X-ray emission could be directly 
caused b y Comptonisat i on of the synchrotron jet compo- 
nent, e.g. iMarkoff et al] (|2005l ). This could explain a num- 
ber of the longer-term corre lation trends in va rious X-ray 
states, e.g., for GX 339-4 (jCoriat et al.ll2009h . What the 
consequence of this would be on shorter time-scales is less 
clear, particularly since the modelling mainly concerned the 



infra-red emission. Fast variabilit y from a jet-like compo - 
nent has been studied in detail by lEikenberrv et al.l (|2008l ) 
for the case of the micro-quasar GRS 1915+105. 

• In any scenario, and indeed in some of the ones above, 
one requires a source of instability to drive emission enhance- 
ments. Such instabilities may naturally arise from the inter- 
actions of magnetic fields and matter, and act as a primary 
channel to emit energy from the system. The presence of 
jets and of high viscosity in the accretion discs are often 
seen as evidence for strong large-scale fields, although more 
direct evidence (such as Zeeman lin es or cyclotron harmon- 
ics) is missing llLovelace et al.ll2009l 'l. Note that for the case 
of SWIFT J1753.5— 0127, humps wer e seen in the spectru m 
that may suggest cyclotron emission (jPurant et al.ll2009l ). 

In summary, there are several plausible ideas for the inter- 
action of low-energy and high-energy emission in accretion 
systems. Typically they have been formulated to explain one 
particular set of observations, and so lack generality or pre- 
dictions. 

As mentioned in some of the models which provide pos- 
sible explanations for our CCFs, multiple spectral compo- 
nent may be involved, some of which span large ranges of 
the available EM spectrum. Both the optical range and 2- 
20keV in X-rays are relatively small windows, and it is not 
at present possible to disentangle spectral components, at 
least for the time-varying portion. In order to separate the 
spectral components within the emission, much wider spec- 
tral coverage is required. Since fast IR devices are now start- 
ing to be available, coupled with wide spectral range X-ray 
satellites such as SWIFT and Suzaku, it is possible to ex- 
tend this work. Naturally, observations of more sources in 
more states will help answer the ubiquity or otherwise of the 
behaviour presented here. 

It is clear that anti-correlations in CCFs are not simply 
a curiosity restricted to a couple of sources. Here we have 
shown their appearance (for at least some of the time) in two 
black hole binaries and two NS binaries. Given the relatively 
sparse sampling and short datasets examined so far, their 
presence here may suggest that this phenomenon is far more 
common than previously thought. Longer coverage during 
various X-ray states is required to test this hypothesis. 
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APPENDIX A: CROSS CORRELATIONS 

In Figurcs rAll to lXMl we show all the CCFs generated in our 
work, in both averaged and dynamic form. There are nine 
panels per figure, with the optical bands red, green, UV 
in columns left-to-right and X-ray bands 2-4 keV, 4-8 keV, 
8-20 keV in rows top to bottom. The exception is SWIFT 
J1753.5— 0127, which was too faint in the UV band for us 
to attempt accurate photometry. Each figure is captioned by 
the Object name, the date, and the RXTE dataset analysed. 

For the average CCF panels, the left-hand y-axis scales 
are in terms of the uncertainty estimate a (see Section 3.2), 
and can be thought of as a rough measure of significance. 
We have kept the y-axis scale the same for all the CCFs of 
a particular object, to enable a fair comparison of changes 
in the CCFs. In each case, the absolute CCF value is also 
given on the right-hand y-axis scale (where the range now 
depends on the noise in each data-set). 

For the dynamic CCF plots, we have kept the same 
colour scaling throughout, since the noise is not necessar- 
ily a constant throughout the observation, and the concept 
of significance becomes somewhat muddied. Note, however, 
that the y-axis scale for the dynamic CCF plots are not the 
same, as this refiects the length in time of each simultaneous 
observation. 




Figure Al. SWIFT J1753.5-0127, 2007-06-12, 93119-02-02-00 
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Figure A2. SWIFT J1753.5-0127, 2007-06-12, 93119-02-02-00 
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Figure A3. SWIFT J1753.5-0127, 2008-08-20, 93105-01-56-00 
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Figure A4. SWIFT J1753.5-0127, 2008-08-20, 93105-01-56-00 
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Figure A5. GX 339-4, 2007-06-13, 93119-01-02 
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Figure A6. GX 339-4, 2007-06-13, 93119-01-02 
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Figure A7. GX 339-4, 2007-06-16, 93119-01-03 
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Figure A8. GX 339-4, 2007-06-16, 93119-01-03 
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Figure A9. GX 339-4, 2007-06-18, 93119-01-04 
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Figure AlO. GX 339-4, 2007-06-18, 93119-01-04 
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Figure All. Sco X-1, 2004-05-17, 90020-01-01-01 
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Figure A12. Sco X-1, 2004-05-17, 90020-01-01-01 
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Figure A13. Sco X-1, 2004-05-17, 90020-01-01-02 
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Figure A14. Sco X-1, 2004-05-17, 90020-01-01-02 
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Figure A15. Sco X-1, 2004-05-17, 90020-01-01-03 
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Figure A16. Sco X-1, 2004-05-17, 90020-01-01-03 
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Figure A17. Sco X-1, 2004-05-18, 90020-01-02-01 
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Figure A18. Sco X-1, 2004-05-18, 90020-01-02-01 
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Figure A19. Sco X-1, 2004-05-18, 90020-01-02-01; immediately succeeds data in Figure A17. 
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Figure A20. Sco X-1, 2004-05-18, 90020-01-02-01; immediately succeeds data in Figure A18. 
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Figure A21. Sco X-1, 2004-05-18, 90020-01-02-02 
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Figure A22. Sco X-1, 2004-05-18, 90020-01-02-02 




Figure A23. Sco X-1, 2004-05-19, 90020-01-03-02 
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Figure A24. Sco X-1, 2004-05-19, 90020-01-03-02 
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Figure A25. Cyg X-2, 2007-10-16, 92038-01-02-00 
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Figure A26. Cyg X-2, 2007-10-16, 92038-01-02-00 




Figure A27. Cyg X-2, 2007-10-19, 92038-01-05-03 
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Figure A28. Cyg X-2, 2007-10-19, 92038-01-05-03 




Figure A29. Cyg X-2, 2007-10-20, 92038-01-06-00 
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Figure A30. Cyg X-2, 2007-10-20, 92038-01-06-00 
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Figure A31. Cyg X-2, 2007-10-20, 92038-01-06-00; follows data in Figure A29. 
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Figure A32. Cyg X-2, 2007-10-20, 92038-01-06-00; follows data in Figure A30. 




Figure A33. Cyg X-2, 2007-10-21, 92038-01-07-01 
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Figure A34. Cyg X-2, 2007-10-21, 92038-01-07-01 
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